Structure-Preserving Transformations of Skew-Hamiltonian/Hamiltonian Matrix Pencils
نویسندگان
چکیده
ityl or positive realness in a VLSI model is an important property Passivity in a VLSI model is an important property to guarantee stato guarantee stable global simulation [3,7]. Existing DS passivity ble global simulation. Most VLSI models are naturally described tests are restrnctive in different aspects. For example, the extended as descriptor systems (DSs) or singular state spaces. Passivity tests linear matrix inequality (LMI) test in [7] has a high complexity for DSs, however, are much less developed compared to their nonof 0(n5) to 0(n6), rendering it prohibitive in testing passivity of singular state space counterparts. For large-scale DSs, the existing high-order DSs, as is usual for VLSI models. The generalized altest based on linear matrix inequality (LMI) is computationally progebraic Riccati equation (GARE) test [8] works only in the limited hibitive. Other system decoupling techniques involve complicated case of admissible (regular, stable and impulse-free) DSs. coding and sometimes ill-conditioned transformations. This paper The contribution of this paper is the formulation of a fast 0(n3) proposes a simple DS passivity test based on the key insight that algorithm for checking passivity of a DS. The key insight is that the sum of a passive system and its adjoint must be impulse-free. when a (possibly impulsive) passive system is added to its adjoint, A sidetrack shows that the proper (non-impulsive) part of a pasthe resulting system, which is again a DS, must be impulse-free. sive DS can be easily decoupled along the test flow. Numerical Numerically efficient and reliable techniques in transforming skewexamples confirm the effectiveness of the proposed DS passivity Hamiltonian/Hamiltonian (SHH) matrix pencils are employed. Aftest over conventional approaches. ter removal of uncontrollable and unobservable impulsive modes, if any, passivity can then be checked through the positive semidefinite of the residue matrix and the proper part of the DS using stan
منابع مشابه
Condensed Forms for Skew-Hamiltonian/Hamiltonian Pencils
Abstract In this paper we consider real or complex skew-Hamiltonian/Hamiltonian pencils λS −H, i.e., pencils where S is a skew-Hamiltonian and H is a Hamiltonian matrix. These pencils occur for example in the theory of continuous time, linear quadratic optimal control problems. We reduce these pencils to canonical and Schur-type forms under structure-preserving transformations, i.e., J-congruen...
متن کاملNumerical Computation of Deflating Subspaces of Skew-Hamiltonian/Hamiltonian Pencils
We discuss the numerical solution of structured generalized eigenvalue problems that arise from linear-quadratic optimal control problems, H∞ optimization, multibody systems, and many other areas of applied mathematics, physics, and chemistry. The classical approach for these problems requires computing invariant and deflating subspaces of matrices and matrix pencils with Hamiltonian and/or ske...
متن کاملComputational Experience with Structure-preserving Hamiltonian Solvers in Complex Spaces
Structure-preserving numerical techniques for computation of eigenvalues and stable deflating subspaces of complex skew-Hamiltonian/Hamiltonian matrix pencils, with applications in control systems analysis and design, are presented. The techniques use specialized algorithms to exploit the structure of such matrix pencils: the skew-Hamiltonian/Hamiltonian Schur form decomposition and the periodi...
متن کاملStructure-Preserving Methods for Computing Eigenpairs of Large Sparse Skew-Hamiltonian/Hamiltonian Pencils
We study large, sparse generalized eigenvalue problems for matrix pencils, where one of the matrices is Hamiltonian and the other skew Hamiltonian. Problems of this form arise in the numerical simulation of elastic deformation of anisotropic materials, in structural mechanics and in the linear-quadratic control problem for partial diierential equations. We develop a structure-preserving skew-Ha...
متن کاملSLICOT Working Note 2013-3 MB04BV A FORTRAN 77 Subroutine to Compute the Eigenvectors Associated to the Purely Imaginary Eigenvalues of Skew-Hamiltonian/Hamiltonian Matrix Pencils
We implement a structure-preserving numerical algorithm for extracting the eigenvectors associated to the purely imaginary eigenvalues of skew-Hamiltonian/Hamiltonian matrix pencils. We compare the new algorithm with the QZ algorithm using random examples with di erent di culty. The results show that the new algorithm is signi cantly faster, more robust, and more accurate, especially for hard e...
متن کامل